Power corrections for jets at hadron colliders 
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We discuss non-perturbative QCD corrections to jet distributions in hadron collisions, 
focussing on hadronisation and underlying event contributions. Using soft gluon re- 
summation and Monte-Carlo modelling we show that hadronisation dominates at small 
values of the jet radius R, behaving as 1/R, while underlying event corrections grow 
with the jet area. This provides a handle to disentangle them and parametrize them 
in terms of measurable QCD parameters, which might enjoy a degree of universality. 

1 Introduction 

With the advent of the LHC particle physics will once again break into new territory at 
the high energy frontier. With 14 TeV available in the center of mass, one would naively 
expect that the dynamics of confinement and low-energy QCD would decouple and, thanks to 
factorization, have a minimal influence on high-pj_ observables (a typical figure of merit being 
A/p± ~ O(10~ 3 )). In general, this expectation is not fulfilled. Even at very high energy, for 
example, no hadronic cross section can be precisely determined without a detailed knowledge 
of parton distributions in the colliding hadrons. Furthermore, a wide range of observables of 
interest for both Standard Model and BSM physics must rely upon a definition of hadronic 
jets and a measurement of the jet energy scale. Experience gained at the Tevatron [5] shows, 
for example, that a 1% uncertainty in the jet energy scale causes a 1 GeV uncertainty in the 
top quark mass determination, and is reflected in a 10% uncertainty in the single inclusive 
jet p± distribution at p± ~ 500 GeVQ. 

Assuming a jet has been reconstructed with an infrared and collinear safe algorithm, 
dependent upon a parameter R defining its size in the rapidity-azimuth plane, the energy of 
the jet will differ from the energy of the hard parton that originated it because of a variety 
of physical phenomena. Radiation from the underlying event and from pileup will spill 
inside the jet cone, increasing its measured energy; on the other hand, radiation produced 
during hadronisation will spill outside the jet cone, leading to a negative correction. It is 
important to realize that, while some of these corrections arc definitely outside the reach of 
perturbative calculations (this is certainly the case for pileup, and to a certain extent for the 
underlying event), hadronisation corrections can be explored with perturbative methods, 
supplemented by soft gluon resummation and power correction technology. These methods 
have been successfully applied to simpler processes such as e + e~ annihilation and DIS (for a 
review, see [4]), where studies of event shapes showed that power corrections to distributions 
can be efficiently parametrized in terms of a limited number of non-perturbative parameters, 
enjoying a remarkable degree of universality 0[6). Here we will apply for the first time these 

a For an early study of the impact of power corrections on jet distributions, see [3]- 
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methods to jet distributions in hadron collisions, a much more challenging environment. We 
consider, as an example, the single inclusive jet p± distribution: we show that perturbative 
methods lead to a prediction for the R dependence of the leading power correction, which 
turns out to be singular (behaving as 1/R) for hadronisation corrections, while, as might be 
expected, it grows as R 2 for the underlying event. We go on to compare the analytic result 
to Monte Carlo models, finding broad agreement as well as some interesting differences in 
the details. 



2 Issues in soft gluon resummation for jet distributions 

Soft gluon resummation for the single-inclusive jet p± distribution was first performed in [7], 
using techniques developed in [8] and [9^1 . Since jet production in hadron collisions 
generically involves at least four colored partons, this is the first case in which nontrivial 
mixing of colour structures occurs. Formally, the structure of the resummed cross section is 



TP ^ 1 

Ej-z — = - exp 
dp j s 



p=l p=J,R 



Tr [HS] . (1) 



The factors £m and £out exponentiate collinear logarithms associated with initial state 
radiation of the colliding partons, and with the measured and the recoil outgoing jets, 
respectively. At the level of power corrections, these factors are expected to generate contri- 
butions of order (A/p±) 2 , associated with jet mass effects. These corrections are negligible, 
and therefore we will concentrate on the contributions of soft gluons emitted at wide angles 
from the jet, embodied in the last factor in Eq. ([1]). Here the trace is taken in the space of 
representations of the color group that can be constructed out of the scattering hard par- 
tons, while H and S are matrices containing hard and soft gluon contributions respectively. 
The exponentiation of soft radiation and the structure of color mixing can be simply under- 
stood [11] by resorting to the eikonal approximation. One can show that all soft logarithms 
can be organized in terms of eikonal colored dipoles, given by 

V tJ (Q, Qo) = f Q a. (*«*>) / d v g r< • Tj , (2) 

where the indices i, j label hard partons, Tj, Tj are the color generators in the corresponding 
representation, and = 2pi- kpj ■ k/(pi -pj) is the transverse momentum with respect to 
the emitting dipole. For n < 4 hard partons, all products of color matrices can be expressed 
in terms of Casimir operators, and thus there is no mixing of color structures. 

Before extending the discussion to power corrections, it should be noted that usage of 
rcsummations such as Eq. |T]) requires great care both in the choice and in the definition 
of the observable. In general jet cross sections, which involve an explicit slicing of phase 
space, are affected by nonglobal logarithms [12] . In the present case, it turns out that the p± 
distribution is a global observable. Nonglobal logarithms will generically spoil Eq. |T]) at NLL 
level, and might influence power correction in a way which is not currently understood. One 
should also note that the choice of jet algorithm may affect the results quite drastically: IR 
safety is a must; furthermore, we will work assuming that the jet momentum is reconstructed 
using four-momentum recombination. 



b A refinement of the implementation of Ref. [7] was recently proposed in [10] 
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3 Radius dependence of power corrections 



In order to analyze the structure of power corrections using the resummation, we consider 
separately each dipole in Eq. |[5J). For the sake of simplicity we place the measured jet 
at zero rapidity, which does not qualitatively affect our results. Let 8^ (k±,rj,<fi) be the 
contribution of a soft gluon with momentum k to the observable, which we normalize as 
£ = 1 — p±/ yS. Note that the contribution is different if the gluon is recombined with the 
jet (<5£ + ), or is left out of it (<5£~). We then construct the shift in the £ distribution due 
to the (ij) dipole by integrating 5^ (k±,r],4>) over the gluon phase space, with a measure 
given by the eikonal dipole, 



.(.ij) 



Pi ■ Pj 

Pi ■ kpj ■ k 



(3) 



Note that 5£ + is integrated inside the jet cone, and 5£~~ outside of it. The result is then a 
function of the jet radius R. 

Different dipolcs give different R dependences, with a transparent physical interpretation. 
The dipole constructed out of the two incoming partons has a leading power correction 
growing like i? 2 , and it is natural to interpret it as the way in which the buildup of the 
underlying event begins to be seen from perturbation theory. Dipoles involving the measured 
jet, on the other hand, behave as 1/R for small R. This behavior arises from gluons which 
are not recombined with the jet: at small R, they are allowed to be emitted very close 
in phase space to the radiating parton, and they begin to see the corresponding collincar 
singularity. All dipole integrals are expressed in terms of a low-energy moment of the strong 
coupling, A(hf), which in principle could be related to the analogous quantity measured in 
e + e~ annihilation. As an example, we consider the contribution of the qq — > qq channel to 
the p± distribution. The dipole comprising the two incoming quarks gives 

A^ in (R) = -^=A(^ F )RJ 1 (R)^-^=A( f i F ) (2R 2 -^R i + 0(R 6 )^ . (4) 

The dipole involving an incoming leg and the measured jet, on the other hand, gives 



A6n-j(*) = ± Afo) (~ - \R + ^R 3 + O (R 5 ) 



(5) 



Notice that, as expected, this corresponds to a negative shift in the p±_ distribution, given 
the definition of £. We remark also that all singular contributions are essentially abelian in 
nature and can be directly collected into an overall shift in the physical distribution, weighed 
by Cf (Ca) for a quark- (gluon-) initiated jet. Terms subleading in R, however, including 
Eq. must be recombined taking color mixing into account, and their contribution is 
more intricate than a simple shift. 



4 Monte-Carlo analysis of hadronisation and underlying event 

A powerful cross-check of the validity of the above discussion is provided by examining 
hadronisation corrections for a range of jet algorithms in Monte Carlo simulations. We 
consider dijet events from both Pythia [13] and Herwig [14] and select events whose highest- 
p± jet at parton-level has 55 < p± < 70 GeV. The selection at parton level is intended 
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Tevatron: 55 < p t < 70 GeV 




hadronisation: f(R) = 1/R 
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Figure 1: Shifts in pt associated with hadronisation and underlying-event contributions for 
the two leading jets in Pythia 6.410 and Herwig 6.510 dijet events. 



to eliminate selection bias associated with non-perturbative effects, as is appropriate for 
comparison with our analytical formulae. For each event we examine the average difference 
in p t for the two leading jets at hadron-levcl as compared to parton level. This difference 
is separated into a hadronisation component and an underlying event contribution, the 
latter 'defined' as that obtained when switching on underlying event (UE) and/or multiple- 
interactions in each Monte Carlo program, using the default parameters of each. The UE 
contribution is expected to be uniform in rapidity and azimuth, scaling as irR 2 for small i?, 
or as 2irRJi (R) for general R with .E-scheme (4- momentum) recombination. 

The upper plots of Fig. [1] show the two non-perturbative components, for two hard 
scattering channels, as a function of R for the kt [E], Cambridge/ Aachen [16], Midpoint^ 
cone [17] and SISCone [18] jet algorithms (all run through Fast Jet [H]). The negative 
result for hadronisation, divergent at small R, is consistent with Eq. ([5]). The rough factor 
of two between the hadronisation corrections for the qq — > qq and gg — > gg is consistent 
with the Ca/Cf ratio expected from Sect. [3] One notes that the cone- type algorithms 
have more negative corrections, differing from the sequential recombination algorithms by 
a term roughly independent of R whose explanation is beyond the scope of the single-gluon 
calculation given above. The UE event contribution to the jet p± is positive, as expected, and 
for the typical range of R studied experimentally, 0.4 < R < 0.7, similar in magnitude to the 

c This algorithm is infrared unsafe and should be thought of as a "legacy" algorithm, shown only for 
historical purposes. 
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hadronisation correction. It is largely identical for all algorithms, and roughly independent 
of the hard scattering channel. 

The analytical R dependence can be studied in more depth by scaling out the expecta- 
tions for hadronisation and UE, as is shown in the lower plots of Fig. [TJ For recombination 
algorithms the result is roughly independent of R, as expected from Eq. ([5]). The normali- 
sation of ~ 0.5 GeV is consistent with the magnitude of hadronisation corrections extracted 
at LEP. In contrast, the normalised UE correction is not constant in R. This can be shown 
to be a consequence of Pythia's implementation of colour-reconncctions between the hard 
partons and the underlying event and when examining models without such reconnections, 
such as Herwig (and also Jimmy [IHj), Fig. [TJi, this effect disappears. One notes that if one 
extracts the UE p± density per unit rapidity (in the figure it is normalised per unit area) 
then it is ~ 6 times larger than the similarly normalised hadronisation correction, suggestive 
of far more violent non-perturbative dynamics. 
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